%线性调频信号LFM
clc;clear;
%D=2*pi*rand;       %相位；rand产生0-1内随机数
f0=0e6;            %中心频率
A=1;               %幅值
T=20e-6;         %采样时间 波形持续时间
B=5e6;             %雷达信号带宽
K=B/T;             %调频斜率
Fs=20e6;Ts=1/Fs;
N=T/Ts;           %采样点数
t=linspace(-T/2,T/2,N);     %时间起始值，中止值，采样点数

%St1=A*exp(1i*(pi*K*t.^2+2*pi*f0*t));
St1=A*exp(1i*(pi*K*t.^2+2*pi*f0*t));
subplot(2,1,1)
plot(t*1e6,real(St1));          %1e6表示10的6次幂；real()St1的实部
xlabel('时间 (us)');
title('Real part of signal');

grid on;axis tight;             %打开网格；将坐标范围设定为被绘制的数据范围
subplot(2,1,2)
freq=linspace(-Fs/2,Fs/2,N);
plot(freq*1e-6,fftshift(abs(fft(St1))));  %fftshift即将零频分量移到频谱中心；abs即取绝对值；fft即做傅里叶变换
xlabel('频率(MHZ)');
title('Magnitude spectrum of signal');
grid on;axis tight;             %打开网格；将坐标范围设定为被绘制的数据范围


